Superfast front propagation in reactive systems with anomalous diffusion 
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We study a reaction diffusion system where we consider a non-gaussian process instead of a standard 
diffusion. If the process increments follow a probability distribution with tails approaching to zero 
faster than a power law, the usual qualitative behaviours of the standard reaction diffusion system, 
i.e., exponential tails for the reacting field and a constant front speed, are recovered. On the contrary 
if the process has power law tails, also the reacting field shows power law tail and the front speed 
increases exponentially with time. The comparison with other reaction-transport systems which 
exhibit anomalous diffusion shows that, not only the presence of anomalous diffusion, but also the 
detailed mechanism, is relevant for the front propagation. 
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The propagation of fronts generated by a reaction- 
transport system has a considerable interest in a large 
number of chemical, biological and physical systems |Q . 
One important model is the advection-reaction-diffusion 
equation (ARD) 



d t e + (u • v)0 = D w 2 e + f(9)/ T , 



(1) 



where Dq is the molecular diffusivity, u is a given incom- 
pressible velocity field and the term f(0)/r describes the 
production process, whose typical time is r. The scalar 
field 6 represents the fractional concentration of the reac- 
tion's products: 9=1 indicates the inert material, 8 = 
the fresh one and < 6 < 1 means that fresh materials 
coexist with products. Usually one considers f(9) with 
an unstable fixed point at 9 = and a stable one at 

Eq. ^) was originally introduced (in the case u = 0) 
by Fisher and Kolmogorov, Petrovskii and Piskunov 
(FKPP) with f(6) =9(1-8) @. Then, the ARD equa- 
tion has been widely studied in the context of combus- 
tion, population growth, aggregation and deposition. 
In absence of stirring (u = 0) it can be shown that a 
front, replacing the unstable state by the stable one, 
moves with a speed vo = 2y/Dof'(0)/T. This result 
is valid whenever f(9) is a positive convex function 
(f"(@) < 0) with two fixed points. For a non-convex 
production term only an upper and lower bound for vq 
can be provided ||. 

A more interesting situation occurs in the presence of 
stirring (u 0): the front propagates with an aver- 
age limiting speed Vf larger than i>o< Under general 
conditions the advection-diffusion equation, i.e., 

dt9 + (u • V)0 = DqSI 2 9, has the same asymptotic fea- 
ture of a diffusion equation. One has that the field (9) 
(obtained with an average on volumes of linear size much 
larger than the typical length of u) at large time obeys 
a Fick equation: 



d t (6) =Y,DtjdU Q ) 



(2) 



The eddy diffusion coefficients, Dfp contain all the (often 
non-trivial) effects of the velocity field u. Eq.(||) states 
that a test particle described by the Langevin equation: 



dx/di = u + \/2D r) 



(3) 



where r\ is a white gaussian noise with (jjj) = and 
{r]i(t)rjj(t')) = 8ijS(t — t'), follows, at large time, a Brow- 
nian motion {(x(t) — x(0)) 2 ) ~ 2D\ x t, where we suppose 
that the first coordinate is the propagation direction. As 
a consequence of the asymptotic diffusive behaviour it is 
possible to show that the front propagates with a finite 
v f <2^/D1J'(0)/T |,|. 

If the conditions for the validity of standard diffusion do 
not hold, anomalous diffusion can be observed, i.e., 



((x(t)-x(0)f)~t 2 



(4) 



with v > 1/2 (super-diffusion) [^Ij- There are at least 
three known mechanisms leading to the anomalous diffu- 

a) an infinite variance of the velocity field u; 

b) an infinite memory, i.e., the velocity- velocity correla- 
tion function has a non integrable tail; 

c) an effective diffusion coefficient increasing with the 
distance between two particles (in the case of relative 
diffusion). 

Mechanism a) is perhaps the simplest one and Levy 
flights belong to this class jl(^,[ll]]. Some deterministic 
maps (e.g., the standard map for specific values of the 
control parameter) can produce super-diffusion according 
to mechanism b) |12|. An important example of mecha- 
nism c) is given by the fully developed turbulence [fl3|| . 
There, the relative dispersion of two particles at distance 
R is described by a diffusion equation with an effective 
diffusion coefficient D e ~ i? 4 / 3 . 
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It is natural to wonder about the effects of the super- 
diffusion for the front propagation. In this letter we dis- 
cuss mechanism a) . The evolution equation of has the 
structure: 



d t = L a 6 + f(6)/T 



(5) 



where L a is a linear operator accounting for the concen- 
tration spreading of test particles evolving according to 
Levy flights with exponent a. An a-Levy flight is an 
independent increment stochastic process, and the dis- 
tribution of each increment is a Levy-stable distribution 



which exhibits a power law tail P a (w) 



, with 



1 < a < 2. The moments behave as (\x(t)\ q ) ~ t q / a for 
q < a and = oo for q > a. So, the role of v in 

Eq. dJ) is played by 1/a (for q < a). Since 1/a > 1/2, 
one can speak of anomalous diffusion. Summarizing, we 
replace the operator — (u • V) + D \7 2 in the Eq. (|l|) by 
L a substituting an a-Levy flight to the Eq. (H) . Only for 
simplicity in the notation and in the numerical computa- 
tion J5|,[l4[ we consider a reacting term which is non zero 
only at discrete time step, when J-form impulses occur: 



f(B,t)= 9(0)5(t-nAt)At 



(6) 



and we introduce the reaction map G{9) = 9 + —g{6) 
governing the evolution of an homogeneous field (i.e., 
without diffusion): 9(t + 0+) = GJ9(t)). The detailed 
shape of G(9) is not important f^JTifl, it is just necessary 
to have an unstable fixed point in 9 = and a stable one 
in 6 = 1. In the following we will present the results for 
the map: 



G{9) = 9/[9 + {l-9) exp(-AVr)] 



(7) 



i.e., the exact solution of the equation d9/dt = 0(1 — 0)/t. 
Noting that, assuming Eq. (JsJ) , between t + + and t+ At 
Eq. (||) reduces to the linear equation <9 t = L a 9, and we 
can determine the field 9(x,t + At) in terms of 9(x,t): 

/ + OC 
AwP a . At {w)9(x -w,t + 0+) 
-OC 

dwP aAt (w)G{d{x-w,t)) (8) 



where P a .At(w) is the probability distribution to have a 
flight of size w in a time interval At. Of course t = nAt, 
where n is an integer number. Let us note that, if one 
assumes the expression (||) for /(0), (§) is an exact re- 
lation and not only an approximation for small At. It 
is important to stress that using (||) one can avoid the 
problem of the precise definition of L a in terms of frac- 
tional derivative |l5|,|l6|. In addition one can study the 
evolution Eq. (||) using a generic distribution P a ,At(w) 
which is in the basin of attraction of the a-Levy stable 



distribution. 

A similar problem has been studied for reaction systems 
driven by a Levy walk p7| . In flLsfl reaction-transport 
processes using non gaussian random walk for the dif- 
fusion process are studied. This approach is similar to 
our one, but the analysis of |1| is always for a class of 
processes which give rise to standard diffusion. In addi- 
tion, in the context of disturbance propagation in chaotic 
extended systems with long-range coupling, Torcini and 
Lepri jl9| studied a discrete space version of (||) with a 
linear shape for G(9). 

Regarding the numerical simulations we study a ID grid 
with open boundary conditions and symmetric initial 
condition with 9{x, 0) ^ in a small region around x — 0. 
We show the results obtained with the reaction map given 
by Eq. (^) and r = 1 . The principal measured observables 
are the quantity of inert material, m(t), and the front 
speed, Vf, defined as: 

f +oc m(t) 
m(t) = / dx9(x,t) v f = lim — ^ . (9) 
J-oo t^ 00 2* 

Simple arguments suggest that if P(w) is steep enough, 
e.g., P(w) ~ exp(— /3|u>|) then the asymptotic behaviour 
of 0(x, t) is the usual one of the standard FKPP equation, 
i.e., 6(x,t) ~ h(x — Vft) where Vf oc y/ (w 2 ) and h(z) ~ 
exp(— jz) for z»l. 




FIG. 1. Front shapes at increasing times in the case of 
P(w) ~ exp(— /3|w|), with (5—1. The inset shows the quan- 
tity of inert material with its asymptotic linear behaviour 
(solid line). We show just one side of the front. 

In this case the inert material increases linearly in time: 
m(t) — const + 2v/t (as for the FKPP equation, see 
Fig. 1). 

More interesting is the case of super-diffusive behaviour 
with a-Levy flight (1 < a < 2). The shape of P(w) 
used in the numerical computation is P(w) = Pq for 
|tu| < Wo an d P(w) = P n \w\~^ 1+a,> for \w\ > w - The 
values of the parameters, P n and Pq, are chosen to guar- 
antee continuity in wq and normalization of P(w). Fig. 2 
shows that at large time develops a power law behaviour 
~ |a;| - ( a+1 ). In addition, the inert material invades the 
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fresh one (6* ~ 0) exponentially fast. This can be seen by 
looking at the total quantity of inert material at time t, 



m(t) 

(see the inset of Fig. 2) . 



(10) 
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FIG. 2. Front shapes at increasing time in the case of 
P(w) ~ |wj~' 1+c *' with a = |. The dashed lines indicate 
the theoretical prediction x~ ( - 1+a \ The inset shows the ex- 
ponential growth of the inert material. The dotted line is 
the asymptotic behaviour proportional to a l ^ 1+a ^ (see the 
main text) and the dashed one is the initial rate which grows 
exponentially as a*. 

The numerical results illustrated above are supported 
by an analytical argument using a linear analysis of the 
tail of 9(x , t) (w hich is expected to be valid for pulled 
dynamics |lJ,E0[) and the theory of the infinitely divisible 
distributions pi| . For 9 around zero, G{9) has a linear 
shape G{9) ~ a9 with a > 1. Plugging this into (||), for 
x 3> 1, one has 



9(x,t) ~ a(P*6)(x,t- 1) 

~ a^P * P * ■ ■ ■ * P * 6)(x,0) 



(11) 



where * indicates the convolution operation. It is well 
known |2lJ] that processes with power law tail P(w) ~ 
|w|~( 1+( ^with 1 < a < 2 are in the basin of attraction 
of the cv-Levy-stable distribution P a (w). Therefore (|lJ 
yields for \x\ ^> 1 and large t 



9{x,t) 



(12) 



in agreement with the behaviour of Figure 2. The growth 
coefficient in Eq. ( |To| ) , c, can be computed with a match- 
ing argument. We expect (and this is confirmed by simu- 
lations) there exists an x value such that 9 ~ 1 for x < x 
and 9 ~ |:e|~( 1+q ) for x > x. We can write m(t) ~ 2i, 
and the value of i can be obtained simply matching the 
0(:M) in (|l|) with 6 ~ 1, i.e., x ~ a*/( 1 + Q ) = exp(^i) 
and therefore c = lna/(l + a), in good agreement with 
numerical results (see the inset of Fig. 2). The transient, 
when max0(cc,O) is small enough, is dominated by the 
reaction term, giving an exponential growth m(t) ~ a*. 



Particularly interesting is the behaviour for P(w) ~ 
w |-(i+q) a > 2. In this case the distribution be- 

longs to the basin of attraction of the gaussian law since 
(w 2 ) < Too (sec ||l]]). In fact, although the probability 
distribution, P(x t ), of the sum of independent random 



variables, Xt = Y^-j 



has a power law tail, the core 



of the distribution behaves as a gaussian, and the tail is 
less and less important as t grows (see Fig. 3). Therefore, 
at first glance, one could expect the same features of the 
standard FKPP equation. 
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FIG. 3. Rescaled probability density functions, P(xt/o-y/t), 
of the sum of independent random variables, Xt = y^*._ t u>j , 

in the case of P(w) ~ |w| _(1+c *' with a = 10 and 
t = 2,4, 8, 16, 32 from top to bottom; a is the standard de- 
viation of P(w). The solid line is the gaussian asymptotic 
behaviour. 

But, in reacting systems the presence of this tail can 
have an important role. In fact for each initial condi- 
tion, 9(x, 0), localized in a small region (e.g., around 
x = 0), already at the first step, the front has a shape 
not steep enough for the usual FKPP propagation (i.e., 
0(x,l) - |:e|-( 1+q )) |0]]. Then, because of the reaction, 
i.e., the instability of 9 = 0, the tail of 9 increases expo- 
nentially in time. As consequence of the gaussian core of 
P{xt) we expect that the bulk of 9 behaves in the FKPP 
way, but, at large time, the exponential growth of the 
tail has the dominant role. 

In Fig. 4 it is shown how, the exponential form of the 
front (which moves with a constant velocity) is overcome 
by a power law tail that grows more and more as the time 
increases, i.e., the initial FKPP behaviour (constant Vf 
and 9(x, t) with exponential decay) is replaced by the ex- 
ponential increasing of the inert material and 9(x,t) with 
power law tail. This result is similar to that obtained 
in p^ ] in the context of the growth of perturbation in 
CML: even a weak long range coupling (e.g, a power law 
with large a) has, at large time, a dramatic effect. 
Summarizing, we have shown how, in a reaction diffu- 
sion system, replacing a standard diffusion with a process 
whose probability distribution has a tail approaching to 
zero faster than a power law, one has the same qualitative 
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behaviour of the usual FKPP system, i.e., exponential 
decay of the front and a finite front speed. 
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FIG. 4. Front shapes at t = 19, 22, 25 • • • where P(w) is the 
same as Figure 3. In the inset it is shown the front speed 
Vf(t) = (m(t + 1) — m(t))/2. The straight lines indicate the 
linear propagation regime, Vf — cost, and the exponential 
propagator regime, Vf ex exp(t lna/(l +a)). 

On the contrary, if the process has the same tail in the 
probability distribution of an a-Levy stable process one 
achieves an exponentially fast propagation of the front, 
instead of the linear one. Moreover the tails of the field 
9 have a power law behaviour instead of the exponential 
decay. As intermediate behaviour we have the case of a 
stochastic process whose increment have a power law tail 
but with finite variance. In this case, initially one has 
the usual FKPP behaviour, but after a while, one has an 
exponential growing of the inert material. 
Let us now compare the above described scenario with 
the case of super-diffusion induced by strong time cor- 
relation of the velocity field. An example is provided 
by the standard map with suitable values of the con- 
trol parameter K fl^| . This advection-diffusion system 
can show super-diffusive behaviour Jl^| in the limit of 
Do — » 0. However in the presence of reaction the front 
speed is finite for each value of Do, because it is bounded 
by U max +VQ, where U ma x is the maximum velocity of the 
test particle, and vq is the front speed in absence of stir- 
ring. Therefore, just the presence of anomalous diffusion 
docs not necessarily imply a non-linear front propaga- 
tion, but also the details of the transport-diffusion and 
reaction mechanisms are important. For example, in the 
system (8) using a reaction map with G(9) — 9 around 
9 = (e.g., ignition) we have always the usual scenario 
shown in Fig. 1 also for 1 < a < 2. 

Such a fact is, in our opinion, relevant for the modeling 
of realistic reaction systems: one has to mimic precisely 
the mechanism that gives rise to the anomalous diffusion, 
beyond the simple power law behaviour of (x 2 (t)). 
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